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SUMMARY 

An  analytical  investigation  into  the  basic  equivalence  of  high-order 
and  reduced-order  aircraft  dynamics  is  presented.  The  need  to  consider 
system  response  to  inputs  represented  by  points  in  the  Laplace  domain  is 
explained  in  terms  of  the  ability  of  a  damped  sinusoid  series  to  mcdel 
general  aperiodic  pilot  inputs.  The  region  of  concern  in  the  Laplace 
domain  is  related  to  pilot  response  time  by  the  time  to  the  initial  peak 
of  a  single  damped  sinusoid  input.  Several  respresentat i ve  high-order 
systems  and  their  lower  order  equivalents  are  examined  for  similarity  of 
time  responses.  A  new  approach  for  equivalent  systems  application  is 
suggested  based  on  matching  time  responses  over  the  Laplace  domain  region. 
Implications  of  the  analysis  to  present  methods  of  applying  equivalent 
systems  and  to  suggested  pilot  compensation  criteria  are  examined. 

The  analysis  leads  to  the  conclusion  that  equivalent  system  para¬ 
meters  are  variables  dependent  on  Laplace  domain  location  and  that  the 
success  of  current  methods  depend  on  that  variation  being  negligible. 

The  further  conclusion  is  reached  that  if  equivalent  system  parameters 
are  allowed  to  vary,  artificial  time  delays  are  unnecessary  and  unde¬ 
sirable  in  achieving  time  response  similarity. 
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INTRODUCTION 


In  the  development  of  the  newest  fighter  aircraft,  increased  stability 
augmentation  and  flight  control  system  sophistication  have  been  employed 
as  a  means  of  maximizing  performance  capabilities.  In  doing  so,  the  form 
of  the  transfer  functions  describing  the  bare  airframe  and  flight  control 
system  have  grown.  In  the  total  system  transfer  function  a  clear  distinc¬ 
tion  between  roots  due  to  feedback  elements  and  those  due  to  bare  airframe 
dynamics  have  often  been  lost.  Since  the  standards  of  MII-F-8785B  are 
based  on  those  root  locations  a  problem  exists  in  how  to  apply  the  require¬ 
ments  and  indeed  if  the  requirements  are  applicable  at  all.  No  extensive 
body  of  data  currently  exists  on  the  acceptability  of  higher-order  sys¬ 
tems  to  pilots.  No  general  consensus  exists  on  proposed  methods  for  deal¬ 
ing  with  systems  of  arbitrary  order.  Thus  both  designer  and  evaluator  are 
left  without  a  clear  cut  means  of  predicting  any  given  system's  acceptab i 1 i ty . 

Recent  attempts  to  deal  with  the  problem  have  included  the  equivalent 
systems  approach  of  reference  (l),  and  the  Neal  and  Smith  criterion  of 
reference  (2).  Although  this  investigation  deals  primarily  with  equival¬ 
ent  systems,  it  will  show  some  important  implications  toward  Neal  and 
Smith's  criterion  as  well. 

Attempts  to  verify  these  criteria  have  concentrated  on  statistical 
sampling  of  pilot  data  and  correlation  of  results  with  the  predictions  of 
each  method.  Since  both  methods  have  been  used  with  some  success,  it  seems 
unlikely  that  statistical  verification  alone  will  resolve  the  differences 
between  them. 

The  following  methods  adopt  a  much  more  analytical  view  defining 
the  concept  of  equivalency  and  developing  conclusions  employing  well  ac¬ 
cepted  linear  control  theory.  The  concept  follows  directly  from  the  Idea 
that  pilots  want  an  airplane  to  fly  like  an  unaugmented  airplane  regardless 
of  its  actual  dynamics.  This  is  the  rationale  under  which  the  existing 
requirements,  or  any  modification  based  on  a  fixed  system  equation,  may 
be  applied  when  the  actual  system  order  may  vary. 
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CONCEPT  OF  AN  EQUIVALENT  SYSTEM 


Whenever  attempting  to  gain  insight  through  analytical  development 
it  is  useful  to  begin  at  the  most  elementary  point.  In  the  case  of  equi¬ 
valent  systems  that  point  is  the  concept  itself  which  may  be  simply  sta¬ 
ted  as  follows:  It  is  possible  that  two  linear  systems  of  different  or¬ 
der  when  stimulated  by  a  common  input  will  produce  highly  similar  outputs. 
If  those  outputs  are  sufficiently  similar  as  to  be  indistinguishable  to 
the  human  operator  of  the  system,  the  operator  will  find  each  system 
equally  acceptable  or  unacceptable.  In  this  event  the  two  systems  may 
be  considered  equivalent,  and  criteria  which  define  the  acceptability  of 
one  may  be  used  to  determine  the  acceptability  of  the  other. 

Current  efforts  in  equivalent  systems  have  taken  the  approach  that 
lower-order  equivalents  may  be  determined  by  finding  a  system  of  prescrib¬ 
ed  form  whose  frequency  response  (juBode)  closely  approximates  that  of  the 
high-order  system  over  a  range  of  frequencies. 

The  form  is  that  of  unaugmented  short  period  and  dutch  roll  for 
longitudinal  and  lateral  directional  dynamics  respectively.  Each  may 
be  modified  by  a  transport  time  delay  to  approximate  high  frequency  lags 
and  obtain  a  more  precise  frequency  response  match.  Since  much  more 
investigation  of  longitudinal  equivalents  has  been  done,  the  analysis 
will  concentrate  on  that  form. 


Longitudinal  Equivalent  System  Form 
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T,  —  ,  Ke,  C  Spe  and  u)nSpe  are  determined  to  minimize  the  magnitude  and 

9  2e 

phase  difference  of  the  high-order  system  ( HOS )  ,  and  low-order  system  (LOS), 
jwBode  between  frequencies  of  0.1  to  10.0  radians.  Computer  programs  have 
been  developed  which  find  the  low-order  equivalent  by  comparing  magnitude 
and  phase  at  a  number  (usually  twenty)  of  frequencies  across  this  range. 
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The  equivalent  system  parameters  are  then  varied  to  minimize  a  squared 
error  function,  termed  a  cost  or  mismatch  function: 

cost  s«2Z^Ma9H0S  "  Ma9Los)2  +  -02  (PhaseH0S  '  PhaseL0S^2 

Magnitude  is  in  decibels  and  Phase  is  in  degrees. 

The  assumption  is  then  made  that  the  resulting  low-order  system  is 
a  true  equivalent  of  the  high-order  system  provided  a  cost  function  of 
10  or  less  is  obtained.  Results  of  early  fixed-base  simulation  experi¬ 
ments  showed  considerable  pilot  dissatisfaction  with  systems  whose  cost 
functions  exceeded  ten  and  therefore  failure  to  obtain  this  level  of  mis¬ 
match  was  itself  considered  reason  for  judging  the  system  unacceptable. 
Time  history  responses  to  step  inputs  were  run  to  confirm  the  assumption 
of  equivalency.  Attempts  to  statistically  correlate  predicted  pilot  ra¬ 
tings  obtained  by  comparison  of  the  low-order  system  to  the  requirements 
of  MIL-F-8785B,  with  actual  pilot  ratings  of  the  high-order  system  have 
shown  good  correlation  in  many  cases  but  unsatisfactory  correlation  in 
some.  The  anomalies  may  be  the  result  of  inconsistent  pilot  ratings, 
but  may  also  indicate  that  the  basic  assumption  linking  frequency  response 
match  with  true  equivalency  does  not  always  hold  for  all  inputs.  The  need 
exists  to  more  clearly  understand  the  link  between  frequency  response 
matching  and  similarity  of  system  time  response.  Since  this  involves 
fundamental  properties  of  linear  systems  the  investigation  is  well  suited 
to  analytical  treatment. 
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EXPANSION  OF  EQUIVALENT  SYSTEMS  INTO  THE  LAPLACE  DOMAIN 

Throughout  this  analysis  several  properties  of  linear  systems  will  be 
employed  which  for  clarity  and  convenience  have  been  substantiated  in  Ap¬ 
pendix  A.  Among  these  is  that  a  transfer  function: 

y(s)  _  k  tt{s-  (am+jbm) } 
xTsT  ~  Tr{s-(an+jbn) } 

when  evaluated  at  a  point  in  the  Laplace  domain,  s  =  a+j-o,  represents  the 
forced  response  portion  of  the  complete  response: 

y(t)  =  GAeat  sinO-ut  +  9)  +  £Bnean1:  sin(bnt  +  'jin) 

■-  ■— ■  n 

forced  response  free  response 

when  the  system  is  subjected  to  the  input: 

x(t)  =  Aeot  sin  wt  ’ 

Since  present  methods  locate  equivalent  system  parameters  based  on 
matching  transfer  function  magnitude,  G,  and  phase  shift,  9,  along  the  ju 
axis,  they  are  enforcing  only  similarity  of  the  forced  response  for  a 
special  class  of  inputs,  undamped  sinusoids.  It  is  not  immediately  clear 
either  that  similarity  of  the  free  response,  or  that  similarity  of  the 
total  response  to  damped  simusoids  are  maintained. 

Of  course,  no  low-order  system  can  approximate  a  higher-order  system 
everywhere  in  the  Laplace  domain.  Forunately  this  is  not  necessary  when  a 
human  pilot  supplies  the  input.  The  principle  is  already  recognized  by 
current  methods  in  that  the  frequency  response  matching  is  done  only  over 
a  limited  frequency  range.  The  upper  limit,  based  on  human  neuromuscular 
lags,  will  later  be  extended  into  the  Laplace  domain.  First,  however,  it 
is  desirable  to  establish  why  consideration  of  the  system  response  to  La¬ 
place  domain  inputs  is  necessary. 

The  basic  concept  of  equivalency  is  based  on  similarity  of  input  and 
output  functions  of  time.  Damped  sinusoids  within  the  pilot's  operating 
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region  are  inputs  to  which  the  system  might  be  subjected.  Therefore  a 
true  equivalent  must  adequately  simulate  the  response  to  such  an  input. 

This  reasoning  alone  is  not  sufficiently  compelling.  It  could  be  argued 
that  since  any  function  possesses  the  Forier  series  expansion: 

x(t)  =  EA  s  i  n  — jj-  t 

that  any  function  may  be  approximated  by  a  sum  of  sinusoids  within  the 
appropriate  frequency  range.  Further,  because  of  the  linearity  of  the 
system,  all  that  is  necessary  is  to  match  the  total  response  to  each  term 
of  the  series  separately  and  the  total  response  will  have  been  matched. 

In  general  pilots  control  aircraft  with  highly  irregular  and  aperiodic 
inputs.  Although  a  Fourier  series  may  represent  a  very  arbitrary  input  sig 
nal,  it  may  do  so  only  over  the  finite  interval:  0  <  t  <  2L.  Outside 
this  interval  the  Fourier  series  repeats;  that  is;  the  function  is  periodi 
regardless  of  the  form  of  the  function  within  the  interval. 

Consider  what  is  meant  by  the  similarity  of  two 'funct ions .  If  two 
functions  of  time  are  identical  then  they  will  have  identical  values  for 
the  function  and  all  its  derivatives  at  any  point  in  time.  If  the 
functions  are  approximate,  only  some  of  these  conditions  will  hold;  that 
is,  the  values  of  the  functions  and/or  their  derivatives  will  match  at 
only  a  finite  number  of  points.  The  number  of  such  similarity  conditions 
which  exist  are  then  a  measure  of  the  accuracy  of  the  approximation.  A 
Taylor's  series  expansion, 

xvt;  i  gpi  „! 

approximates  the  function  by  equating  the  value  of  the  function  and  its 
derivatives  at  t=a.  The  Taylor's  series  representation  is  valid  only  in 
the  vicinity  of  t*a.  To  obtain  an  approximation  valid  over  a  wider  range 
of  time,  the  value  of  the  function  and  its  derivatives  must  be  matched  at 
several  points  over  a  given  time  interval.  Since  it  is  generally  dif¬ 
ficult  to  estimate  higher-order  derivatives  of  a  time  response,  it  is 
most  practical  to  require  that  the  value  of  the  function  and  its  first 
derivative  at  a  point  be  used  as  similarity  conditions.  Obviously 
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guaranteeing  that  only  the  value  of  the  function  matches  the  approximating 
function  at  a  point  says  very  little  about  the  accuracy  of  the  approxi¬ 
mation  at  near-by  points. 

The  general  form  of  the  input  function  which  satisfies  the  differen¬ 
tial  equation  describing  the  system  has  already  been  stated  as: 

x  ( t )  =  Ae at  sin  uit 

and  the  linearity  of  the  system  allows  simultaneous  inputs  of  the  form: 


x(t)  =  -  Aeat  sin  (wt  +  b) 


the  first  derivative  of  which  is: 
dx(t) 


dt 


£Aeat  {jj  cos  ( wt  +  b)  +  a  sin(wt  +  b)  / 


Writing  these  equations  in  terms  of  complex  numbers  simplifies  the  form  to: 
x ( t )  =  eSt  '  -Tp  eS;’;t 


or 


/ .  \  ^  Ae  ^  a 

x(t;  =  £  e51 


and 


AeJ$  . Ae"j^  ....  s+t 


x,(t)  .  i  sest  -  y 


s«e“ 


or 


x'  (t)  »  t 


t 

2j 


where  s  *  a  +  ju>  and  s*  =  a  -  jui. 
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Applying  the  similarity  conditions  at  i  points  yields  the  system  of 
equat ions : 


x(ti)- 

es2ti 

es2 c 1 . .  .  .  es i f 1 

‘  A1e-i<*>1/2 j 

x(t2) 

= 

e^i1; 

es2  '•2 .  .  .  .  es  '  ^.2 

A2e^2/2j 

X (t  i  )  _ 

_  e*i 

es2ti.  .  .  .  eS2l:i  . 

.  Aje-i^i  /2j  _ 

r 

*-» 

X 

"  siesiti  Sje**1!. ..S,es,tl 

A1e-i^1/2j 

X  1  (t,) 

= 

s ^es s.es2t:.  .  .s.esit2 

A2e->  f  2 /2j 

.  x'(t.). 

5  -es  i  i 

L  1  2  1 

A  |  e-i  ®  i  /2  j 

This  system  of  equations  may  either  be  viewed  as  2i  equations  i n*  2 i  complex 
variables;  S;  and  A  { e-J  ^ ' /2  j  ,  or  as  4i  equations  in  the  4i  real  variables; 
oj ,  o»i ,  A;  and  $ ;  .  Either  way  unique  solutions  are  obtainable,  although 
the  solutions  are  not  straight  forward  due  to  the  non  1 i near i t i es  of  the 
equat i ons . 


If  we  consider  only  values  of  s  for  which  s;  =  jwj,  a;  *  0,  we  ac¬ 
tually  eliminate  both  a  and  tu  as  variables  and  the  preceeding  system  of 
equations  becomes  linear  in  the  constants,  Ajej^i/2j. 

The  general  input  form, 

x(t)  *  ZAeat  sin(u>t  +  $) 

with  a  ■  0  becomes: 

x(t)  «  ZA  sin(wt  +  $) . 

Each  individual  term  of  the  series  repeats  its  value  at  tut  *  n2ir.  If 
several  terms  all  repeat  at  the  same  time, 


■  2ir  2  it  2tt 
l. —  =  nt —  »n  —  - 


<*>1 


U)2 


'■U  3 


etc. 


I 
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which  implies, 


uii 


m 
1  ’ 


ui2 


_n 

m’ 


etc . 


This  implies  that  the  ratios  of  all  frequencies  in  the  series  are  express- 
able  as  ratios  of  integers;  i.e.,  rational  numbers.  Since  any  frequency 
may  be  approximated  by  a  finite  decimal,  the  periodicity  of  the  series 
holds  to  any  desired  accuracy  for  any  combination  of  frequencies.  If 
the  series  repeats  over  the  interval  2L,  all  frequencies  are  integral 
multiples  of  a  single  base  frequency. 


.  TT  Tt  IT 

idi  ■  l-j-,  cu2  -  fr>j— ,  J3  =  rrj-  -  -  -  etc. 

Thus,  with  a  =  0  the  general  input  form  becomes  that  of  the  general  Forier 
series , 

x(t)  =  £  A  sin  (-^t  +  4>)  , 

in  order  to  avoid  the  Fourier  series  repeating  within  the  time  interval  of 
interest, 


8ey ond  this  the  choice  of  L  is  arbitrary,  but  once  chosen  sets  the  values 
of: 

s  ”  1  <  n  <  i 

Therefore,  Sj  having  been  determined,  half  the  variables  in  the  matrix 
equations  are  eliminated,  and  only  half  of  the  2i  similarity  equations 
may  be  satisfied.  The  damped  sinusoid  series  then  represents  a  better 

approximation  than  does  the  Fourier  series  for  a  finite  number  of  terms  be¬ 
cause  it  can  satisfy  twice  as  many  similarity  conditions.  These  reasons 
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are  considered  sufficient  to  justify  consideration  of  system  response 
to  inputs  in  the  Laplace  domain.  It  remains  to  more  clearly  define  the 
region  of  the  Laplace  domain  within  which  human  pilots  are  capable  of 
operat i ng . 

It  is  generally  conceded  that  a  pilot's  reaction  capabilities  pro¬ 
vide  an  upper  limit  to  the  frequency  range  considered  for  j ^  Bode  match¬ 
ing  techniques.  The  problem  exists  of  how  to  extend  this  principle  into 
the  s  plane  in  order  to  consider  system  response  to  aperiodic  inputs. 

The  frequency  of  a  sinusoidal  input  is  easily  relatable  to  any  fraction 
of  a  cycle  deemed  to  be  most  significant.  The  only  general  similarity 
between  an  undamped  and  highly  damped  sinusoidal  is  that  both  possess  an 
initial  maximum  or  peak. 


With  the  interpretation  that  the  initial  peak  represents  the  time  to 
initiate  removal  of  a  previous  control  motion,  it  may  be  related  to 
response  time.  The  initial  peak  then  represents  a  convenient  means  of 
extending  the  upper  frequency  limit  over  the  s  plane.  The  general  in¬ 
put  form  is: 

x(t)  *  AeCTt  sin  jt 

At  the  first  peak,  its  first  derivative  vanishes: 


dx(t) 

dt 


which  reduces  to: 


or : 


Ae 


at 


{a  sin  +  w  cos  ^t  } 
o  o 


a  sin  wt  +  u>  cos  wt 
o  o 


t  ■  —  tan 

O  a) 


-i 


pi) 


Note  that  at  a  *  0  the  time  to  the  peak  becomes  the  quarter  cycle  time; 

'o-fc 
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which  for  the  customary  upper 

t  = 

o 

The  pilot  operating  region  is 

a 

Note  that  as  u  +  0,  tan  wtQ  -*■ 

1  im 

■w  -*■  0  0  = 


frequency  limit  of  10  radians  per  second 


20 

then 

s  i  n 


=  .157  sec. 
defined  by  the 

tan  Ut0) 

.utg  iwtp  and : 


con  tou  r : 


~  =  -6.366  sec~ ] 
to 


IHECiijlivG  f„QE 


i>c  1  Fl 
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TIME  RESPONSE  MATCH  I  MG 

The  final  test  of  any  equivalent  is  its  ability  to  approximate  the 
total  time  response  to  an  arbitrary  aperiodic  pilot  input.  This  is  equi¬ 
valent  to  matching  time  responses  over  the  defined  region  in  the  s  plane. 
There  is  therefore  a  direct  implication  that  equivalent  system  parameters 
should  be  determined  by  maximizing  the  similarity  of  the  total  time  re¬ 
sponses  to  such  inputs.  We  may  consider  applying  to  the  output  the  same 
similarity  conditions  we  have  already  applied  to  the  input.  The  general 
output  form  is: 

y(t)  =  GAeat  sin(wt  +  9)  +  E8neant  sin  (bnt  +  >i»n) 
or  in  terms  of  complex  variables: 

y(t)  =  GAe-l®est  -  GAe"J®es  t  +  £Bnej'l'n ernt 
and  its  derivative: 

y  1  ( t )  *  sGAe-i®est  -  s  *GAe”-i®es"t  +  Z  rnBne-i^nernt 


where  s  ■  o  +  jw,  s* 


rn  =  an  ±  jbn,  and  where: 


r  j q  _  K  tt{s  -  (am+j  bm) } 
ir(s  -  (an+jbn) } 


at  s  =  o  +  jw 


The  number  of  variables  in  these  equations  is  directly  dependent  on  the 
order  of  the  transfer  function  and  is:  Number  of  variables  =  1  +  m  +  2n. 
It  is  possible  therefore  to  enforce  this  number  of  similarity  conditions 
on  the  low-order  system  output.  It  is  also  possible  to  consider  other 
means,  such  as  a  minimum  squared  error  method,  of  enforcing  time  response 
similarity.  In  addition,  considerable  judgment  may  be  involved  in  the 
choice  of  times  at  which  to  enforce  those  conditions.  These  difficulties 
will  not  be  critical  to  the  most  significant  conclusion  of  the  analysis* 
but  for  the  present  will  be  dealt  with  by  making  certain  assumptions 
about  the  pilot's  perception  of  the  output. 

If  it  is  assumed  that  human  pilots  are  not  sensitive  to  small  dif¬ 
ferences  in  the  system  response  but  only  to  certain  gross  characteris¬ 
tics,  we  may  set  similarity  conditions  which  duplicate  them.  Based  on 
pilot  comments  on  high-order  systems  contained  in  references  3  and  k, 

*  see  page  55 
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the  most  common  judgments  made  on  the  systems  seem  to  fall  into  three 
catagor i es  : 

1:  Response  magnitude  or  "sens i t i v i ty" 

2:  Response  lags  or  sluggishness 

3:  Annoying  residual  oscillations  or  tendency  to  PIO. 

For  well-damped  inputs  such  as  that  in  Figure  1-a,  the  pilot  will  most 
likely  interpret  the  first  peak  of  the  response  as  the  system  response 
to  his  input  and  the  long  term  response  as  characteristic  of  the  re¬ 
sidual  oscillations.  Matching  the  timing  and  magnitude  of  the  initial 
peak  would  guarantee,  under  these  assumptions,  similarity  of  system  sen¬ 
sitivity  and  speed  of  reaction  as  perceived  by  the  pilot.  Since  the  pilot 
could  consider  the  residual  oscillations  an  annoyance,  it  is  probably  most 
significant  to  match  the  time  at  which  they  damp  out. 

The  first  "n"  similarity  conditions  are  always  set  by  the  require¬ 
ment  that  the  system  be  initially  undisturbed.  The  values  of  the  am¬ 
plitudes  and  phasing  of  the  free  response  modes: 

BneJ’('nernt  rn  *  an  ±  jbn 

are  set  by  requiring  that  the  value  of  the  output  function  and  its  first 
n-1  derivatives  be  zero  at  time  zero.  Considering  the  first  over  second 
order  equivalent  system  form  currently  in  use,  these  conditions  are: 

GAe J 9  -  GAe-J9  +  B1eJ1('i  +  B^-^2  =  0 
sGAe-i9  -  s*GAe'J 9  +  r1B1eJ’l'i  +  r2B2eJ4'2  =  0 

Applying  similarity  conditions  at  the  peak  time,  t^,  and  the  subsidence 
time,  t2, 

GAeJaestl  -  GAe’-^e5*1^  +  B1e*1eriti  +  B2e*2er2tl  =  Ypk 
sGAe-*8estl  -  s*GAe'-i9es*ti  +  r1B1e'*'1eritl  +  r2B2e'i'2er2t2  =  0 
GAe-i9est2-  GAe’j9es*t2  +  B^e^2  +  Bje^e^2  =  Ym,n 
sGAeJ0est2  -  s*GAe'J9es*t2  +  r1B1e'^1erit2  +  r2B2e’l'2er2t2  =  Y'min 
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where : 

A,  s,  and  s*  are  known  because  :~.e  incut  is  known 
Y  r  ’  anc*  c'  are  determined  cv  the  hign-order  system  response 
which  is  known 

Ym;n  and  Y'm;n  are  values  of  the  function  ana  its  slope  which 
define  the  subsidence  time  t  .  They  are  small  but  not  exactly  zero. 
G,  f ,  r .  ,  r,,  B  ,  e  "  * ,  and  S.e":  are  variables. 

The  similarity  equations  may  be  written  more  conveniently  in  Matrix  form  as: 


0 

A  A  11 

GeJ  9 

0 

As  As*  r  r 

1  2 

Ge"J  9 

Ypk 

Aest  l  Aes';ti  eriti  er-ti 

B1eJ'a'i 

0 

sAesti  s*Aesti  r1eriti  r2er2,:i 

B2eJ^ 

Ymin 

Aest2  Aes':t2  erit2  er2t2 

Y'min 

sAesta  s*Aes"ta  rTerita  r,er2ta 

Although  highly  nonlinear,  this  system  of  equations  has  a  solution  best 
obtainable  by  numerical  methods.  Having  obtained  a  solution  for  the  trans¬ 
fer  function  magnitude,  phase,  and  natural  response  modes  (i.e.,  pole 
locations),  the  numerator  zero  location  and  gain  are  determined  by: 

K(s  -  a 3)  =  GeJ'9(s  -  rj  (s  -  r2) 

The  high  order-system  could  be  approximated  by  lower  than  a  first 
over  second-order  system  but  only  if  some  of  the  similarity  conditions  are 
sacrificed.  Higher-order  equivalents  would  produce  better  time  response 
matches  but  not  different  pilot  ratings  under  the  assumptions  just  made. 

Figure  1-b  represents  a  hypothetical  high-order  system  response  to 
a  highly-damped  input,  Figure  1-a.  Figures  1-c  to  1-h  represent  possible 
low-order  system  responses  to  successively  higher-order  equivalents.  These 
show  how  the  equivalents,  by  enforcing  progressively  more  similarity 
conditions,  improve  the  time  response  matches.  For  instance,  Figure  1-c 
demonstrates  that  a  pure  gain,  having  only  one  variable,  is  capable  only 
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of  matching  the  peak  magnitude  of  the  high-order  system  response.  Figure 
i-g  demonstrates  that  the  first  over  second-order  equivalent  is  the  sim¬ 
plest  form  capable  of  imposing  similarity  of  the  function  initially,  at 
the  first  peak,  and  at  response  subsidence  simultaneously. 

When  the  time  delay  is  added,  as  in  Figure  1  —  h ,  a  time  scale  shift 
is  effected  and  the  similarity  conditions  which  were  applied  at  t  =  0, 
t  =  tx,  and  t  =  t2  are  now  applied  at  t‘  =0,  t[  =  tx  -  T,  and 

=  t,  -  T.  The  result  is  that  although  the  variable  T  was  acquired, 
the  initial  conditions  are  applied  at  t  =  T  rather  than  t  =  0  and  those 
two  similarity  conditions  are  lost.  The  addition  of  the  time  delay  re¬ 
sults  in  a  net  loss  of  similarity  conditions  and  detracts  from  the  overall 
.similarity  of  the  time  response.  If  it  is  sufficiently  small,  the  delay 
may  not  significantly  detract  from  the  initial  response  similarity  but  will 
in  general  change  the  other  equivalent  systems  parameters. 
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SAMPLE  HIGH  ANO  LOW  ORDER  SYSTEMS 


To  this  point  in  the  analysis  concepts  and  principles  have  been  dis¬ 
cussed  which  will  serve  as  a  framework  for  the  analysis  of  specific  high- 
order  systems  and  their  equivalents.  The  high-order  systems  were  selec¬ 
ted  from  among  those  of  the  LAHOS  experiments  of  reference  2,  while  their 
lower  order  equivalents  were  obtained  from  reference  1.  These  systems 
were  the  LAHOS  1-4,  LAHOS  1-C,  and  LAHOS  6-2  configurations.  LAHOS  1-C 
was  considered  typical  of  those  in  the  references,  LAHOS  1-4  represents 
an  inadequate  form  according  to  the  previous  development,  and  LAHOS  6-2 
represents  the  highest  order  reduction  in  the  references.  All  possess 
good  Bode  matches  with  cost  functions  less  than  10.  The  notation  used 
will  be : 

(TS+1)  =  [T]  ;  (S2/'4+2'/unS+l)  =  Ln>  cl 

or: 

( S  -  c )  =  (a);  (S‘  -  2jS  +  a2  +  mZ)  =  (c , 

CONFIGURATION:  LAHOS  1-4:  0/2nd  Order  Equivalent,  Cost  Function  9.0: 

High  Order  System: 

j_  _  [-1-4] _ _ 

Fs  ~  [26.6,  0.6K75.0,  0.7JL1.0,  0.74)1-0.5] 

or : 

9  1 .0647  x  107  (-0.714) 

Fs  "  (-15.6,  20.8) (-52.5,  53 . 56)  (-0 . 74 ,  0.6726) (-2.0) 

Low  Order  System: 

0  1.08 

Fs  =  [1.56,  0.74J 

or : 

_9_  2.628 

Fs  "  (-1.15,  1 .05) 
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Time  Delay 


-0.06  seconds 
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CONFIGURATION:  LAHOS  1-C:  L3  Free,  Cost  Function  8.9: 

High  Order  System: 


±  _  [-1.4][-0.2] _ 

Fs  [26.0,  0.6] (75.0,  0.7JL1.0,  0.74JI-0.1] 


or: 


J_  =  1.0647  x  107  (-0.714) (-5-0) _ 

F s  '  (-15.6,  20.8) (-52.5,  53 . 56) (-0. 74 ,  0.6726) (-10.0) 

Low  Order  System: 

i.  _  0. 99 [-2 . 63] 

Fs  [0.86,  1.04] 
or : 

0  _  1-9257  (-0-38) 

Fs  (-1.14) (-0.6487) 

Time  Delay  =  -0.037  seconds 


CONFIGURATION:  LAHOS  6-2:  L^  Fixed,  Cost  Function  0.97 

High  Order  System: 

J_  _  [0.5] [0.43] [0.06] [1 .4] _ 

F s  ~  [0 .2] [0 . 1 ] [ 1 . 1 ] [ 1 . 9,  O.65] [26.0 ,  0.6]  [75.0,  0.7] 

j_  _  1.1269  x  IQ7  (-2.0) (-2.33) (-16.67) (-0-71) _ 

F s  "(-5.0T C-10.0J (-0.91  ) (-1 .235,  1.444) (-15.6,  20.8) (-52.5,  53-56) 


Low  Order  System: 


1-  -97  [1.4] 

Fs  [1.74,  0.78] 

9_  4.111(-0.714) 

Fs  =  (-1.357,  1.089) 


Time  Delay  =  -0.084  seconds 
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DISSIMILARITY  OF  THE  FORCED  RESPONSE 

Since  current  equivalent  systems  methods  match  the  forced  response, 
transfer  function  magnitude  and  phase,  along  the  jw  axis,  we  might  consider 
whether  this  principle  is  extendable  into  the  S  plane.  If  similarity  of 
the  forced  response  per  se  is  a  necessary  condition  to  similarity  of  the 
total  response,  then  matching  the  forced  response  is  a  general  principle 
extendable  over  the  appropriate  S  plane  region.  Since  equivalent  systems 
parameters  are  currently  being  determined  to  minimize  the  cost  function 
along  '  '  jw  axis,  it  is  reasonable  to  believe  that  no  such  similarity 
exists  away  from  the  jw  axis.  Tables  one  through  six  compare  high  and  low 
order  transfer  function  magnitude  and  phase  angle  of  the  three  LAHOS  con¬ 
figurations  under  consideration  over  a  frequency  range  zero  to  ten  radians 
and  sigma  values  from  zero  to  -8.5- 

The  phase  angles  in  these  tables  occasionally  appear  positive.  The 
phase  angle  is  calculated  by  an  inverse  tangent  function  which  returns 
angles  between  ±180°.  Positive  phase  lags  then  actually  represent  the 

an9,e:  9  =  -360°  =  9 comp  9comp  >  0 

For  the  low-order  equivalents,  values  greater  than  -180°  occasionally  appear 
because  the  phase  shift  due  to  the  time  delay  is  added  to  the  inverse  tan¬ 
gent  function. 

From  these  tables  it  is  apparent  that  no  similarity  exists  anywhere 
but  along  the  line  a  *  0.  Therefore  if  similarity  of  the  forced  response 
is  a  valid  general  principle  there  would  be  considerable  differences  in  the 
total  time  responses  to  damped  sinusoidal  inputs.  Figures  2,  3,  and  4 
have  been  plotted  to  more  clearly  show  the  forced  response  variation  of  the 
three  configurations  at  one  radian  frequency  across  this  damping  range.  The 
sigma  limits  dictated  by  extension  of  pilot  reaction  times  into  the  S  plane 
are  also  shown. 
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CONFIGURATION:  LHHOS  1-4  HIGH  OPDEF  SYSTEM 


TRBNSFEP  FUNCT I ON  6A I N 


GMA  ?  OMEGA 

.  0  0 1  0 

1 . 0  0  0  0 

3. 0000 

5. 0000 

8. 0000 

1  0 .  0  O  i'i 

.  0  0  0  0 

.  9996 

1 . 0401 

•  cl  6  d  \ 

.  1  056 

.  04  35 

.  0384 

-.50 0 0 

.  S  08  0 

d ■ o9c r 

m  J  U  ll  4 

.  1138 

.  0458 

.  0388 

-1.0 0  0  0 

1 . 6484 

3. 1 051 

.  3316 

.  1807 

.  0430 

.  03 1 1 

-1.5000 

4. 7048 

3.3081 

.3441 

.  1853 

.  0488 

.  0334 

-£.  0  0  0  0 ♦ C 

5. 4364 

1 . 8033 

.  3388 

.  1874 

.  0514 

.  0335 

-8 .  5  0  0  0 

3. 30c 4 

1  • 

.3181 

.  1878 

.  0585 

.  0345 

-3. 0000 

1 . 38 1  6 

.  8308 

.  8806 

.  1848 

.  0533 

.  0354 

-3.5  0  0  0 

.  3  03 1 

.  6383 

S 

.  1807 

.  0537 

.  0361 

-4 . u  0 0  0 

.5331 

.4605 

“r  O  ”7  CT 

.  1153 

.  0537 

.  0367 

-4 . 5  0  0  0 

.  3844 

.3471 

.  1  838 

.  1  088 

.  0535 

.  0371 

-5. 0000 

.  £838 

.8716 

.  1 738 

.  1  087 

.  0588 

.  0374 

-5.50 0  0 

.3381 

.  8181 

.  1585 

.  086 1 

.  0581 

.  0375 

-b. 0000 

.  1885 

.1311 

.  1344 

•  0397 

.  05 1 1 

.  0375 

-6.5  0  0  0 

.  1534 

.  1587 

.1181 

.  03  2:3 

.  0500 

.  0374 

-7. 0000 

.  135  0 

.131  0 

.  1  068 

.  0778 

.  0438 

.  0378 

-7.5  0  0  0 

.1167 

.1138 

.  0853 

.  0736 

.  0474 

■  Li  *:  6 

-3. 0000 

.  1  03  3 

.  1  0  0 1 

.  0860 

.  0677 

.  0460 

-  0366 

-8.5 0  0  0 

.  0806 

.  088  0 

.  0780 

.  0633 

.  0446 

.  0361 

TRANSFER  FUNCTION  PHASE  LAG 


1 5  OMEGA 

.  0  0 1  0 

1 .  0  0  0  0  3 .  0  0  0  0 

5. 0000 

8. 0000 

1  0 .  0  0  0  0 

.  0  0 

-.  04 

-65.81 

-141.87 

-163. 08 

1 63 . 8  0 

1 56 . 58 

-.50 

.  17 

-85. 05 

-158. 01 

-178.87 

161. 88 

151.0 0 

-1. 00 

1  r‘  9  >  y  u 

- 1 68 . 88 

-177. 58 

168.83 

1 54 . 86 

145.45 

-1.50 

178.88 

148. 87 

163.48 

1 58 . 6 1 

1 47 . 88 

1 38.88 

-8.  00 

80.  03 

1 16.  1 1 

144. 88 

147.38 

141 . 07 

134. 43 

—3 . 5  0 

.  14 

84 .  78 

1  d  r‘  •  9 

136. 48 

1 34 . 85 

188. 01 

-3.  00 

.  08 

68 . 6 1 

118.51 

1 £6 . 08 

1 £7 . 58 

183.  7  0 

-3.  50 

.  05 

48 .  46 

88.  07 

1 16.30 

181.18 

118.58 

-4.  00 

.  04 

38 .  1  £ 

37.54 

107.33 

114.81 

113. 50 

-4.50 

.  03 

38.61 

f'  •  T’  5 

88 . 8 1 

1 03 . 87 

1 09 • 66 

-5 .  0  0 

.  03 

£7 . 84 

68.  47 

81 . 34 

1 03 . 33 

1  04.  01 

-5.50 

.  03 

34.31 

6£ .  47 

84.51 

88.  01 

88.  58 

-6.  00 

.  03 

£  1 . 37 

56.  53 

78 . 36 

83 .  03 

85 . 4  0 

-6.  50 

.  OS 

1 8 .  08 

51.43 

78.85 

88 .  38 

81.45 

-7.  00 

.  03 

1 7 .  ££ 

47.15 

67. 8£ 

84 .  06 

S  t'  •  75 

-7.50 

.  03 

15.67 

4  3 .  44 

63.  58 

30.  07 

84. 30 

-6 .  0  0 

.  01 

14.37 

4  0.83 

58 . 6  0 

76.41 

81.1  0 

-8. 50 

.01 

13. £7 

37 . 44 

56.10 

73 .  05 

78.  1 7 

TABLE  1  - 

TRANSFER 

FUNCTION 
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CDN^  I  S'JR’RT  I  DM :  LRHO.-  1-4  0 /SND  ORDER  ECUIVRLENT 


TR'RNiFEP  FUNCTION  GRIN 


S I G M h ! OMEGA 

.  001  0 

1.0  0  0  o 

3. 0000 

5 .  0  0  0  0 

S. 0000 

1 0. 0000 

•  0000 

1 . 083? 

.3713 

.  8757 

.  1  037 

.  04  03 

.  0868 

-.50 0 0 

1 . 7833 

1 . 8745 

.3117 

.  1  073 

.  04 1 5 

.  0865 

-1.  0  0  0  0 

8. 33*u 

3. 0868 

•  3  3 1 6 

.  1  033 

.0413 

.  0866 

-1.5000 

8. 1453 

3.5748 

.  3863 

.  1  034 

.  0417 

.  0865 

-8. 0000 

1 . 4400 

1 . 3308 

.  8385 

.  1  065 

.  041 3 

.  0864 

-  2. .  5  0  0  0 

.  8335 

.  8536 

.1016 

.  04  06 

.  0861 

-3. 0000 

.  5308 

.5143 

.  8 1 36 

.  0358 

.  0396 

.  0857 

—  3.5  0  0  0 

•  3967 

.  3585 

.  13  33 

.  083 1 

.  0  383 

.  0851 

-4. 0000 

.  8843 

.  8686 

.  15  37 

.  OS  07 

•  0369 

.  0845 

— 4 .  50  0  ij 

.  81 38 

.  1337 

.  1830 

.  0734 

.  0353 

.  0833 

-5. 0000 

.  165  0 

.  1 565 

.  1  030 

.  0664 

•  0336 

.  0831 

-5. 5000 

.  1318 

.  1 856 

.  0388 

.  060  0 

.  0313 

■  IJilklu. 

-6. 0000 

.  1  067 

.  1  083 

.  0736 

.  0548 

.  0308 

.  08 1 4 

-6.50  0  0 

.  0334 

.  0857 

.  06  38 

.  0483 

.  0885 

.  08  05 

-7. 0000 

.  0744 

.  0785 

.  0533 

.  0448 

.  0868 

.  0137 

-7.5  0 0 0 

.  0634 

.  068  0 

.  0585 

.  0401 

.  0853 

.  0 1 88 

0000 

.  0547 

.  0537 

.  0464 

.  0364 

•  iJ  c!  r 

.  0173 

-8.50  0  0 

.  0477 

.  0463 

.  0418 

.  03  31 

.  0883 

.0171 

TRANSFER  FUNCTION  PHRS 

sigmr; omegr  .0010 

E  LRG 

1 .  0  0  0  0 

3. 0000 

5. 0000 

3. 0000 

1  0 .  0  0  0  0 

.  00 

-.  06 

—6 1 . 66 

-143. 33 

-170.13 

-130. 86 

-801 . 1 1 

-.50 

-.  05 

-71.45 

-168. 76 

-181.71 

-138. 05 

-8  06 . 86 

-1.0  0 

-.  08 

— 70 . 88 

-133.73 

-133.53 

-805. 38 

-818. 64 

-1.50 

.  03 

63 . 74 

154. 57 

1 54 . 44 

147.40 

141.57 

-8.  00 

.  05 

6  0 . 68 

134.83 

1 48 . 67 

140.  16 

1 35 . 3  0 

-8.50 

.  05 

51 . 07 

1 16.56 

131.36 

1 33 . 08 

130. 08 

-3.  00 

.  04 

48 . 35 

101.64 

180.71 

1 86 . 04 

184. 44 

-3.  50 

.  04 

36 . 44 

83.  85 

no.  S3 

113. 86 

113. 30 

-4.  00 

.  03 

31.88 

78.  33 

101.73 

118. 78 

113.50 

-4.50 

.  03 

8  r‘  .  1  r' 

70.83 

33 . 54 

1 06 . 45 

1 08.84 

-5.  00 

.  08 

83.85 

63.  00 

86.  07 

1 00, 47 

103.  14 

-5.50 

.  08 

81. 14 

56.73 

73.34 

34.  78 

38.81 

-6 .  0  0 

.  08 

18.  88 

51.45 

73.85 

83.40 

33.47 

-6 . 5  0 

.  08 

16.33 

46.84 

67.  76 

84 .  38 

88.  38 

-7.  0  0 

.  08 

1 5 . 38 

48. 88 

68.80 

73.53 

84 . 56 

-7.50 

.  01 

14.  00 

33.83 

53.31 

75.  03 

80.  33 

-8.  00 

.  01 

18.80 

36.  17 

54.83 

70.  73 

76.40 

-8.50 

.  01 

11.76 

33 . 4  0 

50.58 

66.  81 

78.60 

NADC  ~  792  3 1 -60 


CONF I GUR  AT I DN :  LhHQ:  1-C  HIGH  ORDER  SYSTEM 

TRANSFER  FUNCTION  GRIN 


isma; omega 

.  0  0 1  0 

1.  00  0  0 

v  3. 0000 

5. 0000 

8. 0000 

1  0 .  Ci  i'i  i'i 

.  UUUO 

.  999*6 

1 .  1800 

.  5290 

.  3597 

.  2644 

- .  5  u  u  u 

.  5748 

2. 0816 

.  55  06 

•  i"  c.  r' 

.2761 

.  2396 

-1.0  i.i  i.i  o 

.  7500 

1 . 9994 

.5527 

■  Jc-c!  i 

•  £87  6 

.  25  06 

-1.5  ij  U  0 

>  9696 

1 . 0975 

.5353 

-  3895 

.  2990 

.  2618 

-2 .  t.i  o  o  o 

.  7521 

.  7461 

.5047 

.  3938 

.3104 

.  £733 

-2.50  0  Cl 

.  55  04 

.5504 

.  4692 

.  3964 

.  3218 

.  2850 

-i. 0000 

•  3*976 

.4163 

.  435 1 

.  3984 

>  -_r*  s! .  -! 

.  2971 

-  ji .  5  0  ij  ij 

.  £?80 

.3159 

.  4067 

.401  1 

.  3450 

.  3  094 

-4. 0000 

.  1  7  77 

•  C-.-94 

■  o  8  6  6 

.4055 

.  3572 

.  3222 

-4.5  i.i  ij  0 

•  Uc'  r‘  4 

«  1 8  6  9 

.  3770 

.4127 

•  3699 

.3353 

-5.  0  i.i  0  0 

.  0002 

.  1 684 

.  3795 

.  4234 

■  -•  •-< 

.  3488 

-5.5000 

.  0902 

.  1934 

.  3952 

.  4383 

.  3975 

.  3627 

-b.  ij  ij  i.i  i.i 

.  1895 

.2561 

.4249 

.  4576 

.4123 

.  377  0 

-6 . 5  ij  ij  Ci 

.  3055 

•  3487 

.  4687 

.4812 

.  4279 

.  39 1 7 

-7. OOUO 

.  4499 

.  4?££ 

.  5264 

.  5  089 

.4440 

■  4  066 

-7 . 50  0  0 

■  64c  1 

.  6365 

.5971 

.  5397 

■  46  04 

.  42 1 6 

-8 .  0  ij  iJ  0 

.  9205 

.8615 

.  6790 

.  5727 

.4770 

.  4367 

— d .  5  Ci  ij  i.i 

1 . 3737 

1 . 1813 

.  7679 

.  6062 

.  4934 

.4516 

TRANSFER  FUNCTION  PHASE  LAG 


SIGMA: OMEGA 

.  001  0 

i  1.0  0  0 

0  3 .  0  0  0  0  5 .  0  0  0  0  8 .  0  0  0  0  1  0 .  0  0  0  0 

.  ci  ij 

- .  o  0 

—33 . 64 

-71 . 29 

-81.39  -95.90  -106.35 

-.50 

.  22 

-54.84 

-79. 41 

-1.0  0 

179.86 

-117.30 

*br • 5 W 

-89.09  -100.36  -110.07 

-1 . 5  i.i 

-179.98 

—  1  37 . 36- 

-94. 88 

-92.56  -102.48  -111.87 

-2 .  0  ij 

-179.97 

-142. 58 

-1 00. 56 

-95.58  -104.49  -113.61 

-8.50 

-179. 96 

-144.51 

-104.27' 

-98 . 06  - 1 06 .38  -115. 29 

-3.  00 

-179.96 

-143.95 

-105.94 

-99.95  -108.14  -116.91 

—•3. 5  0 

—1 79. 96 

-140.28 

-1 05. 70 

-101.27  -109.78  -118.46 

-4.  00 

-179.94 

-131. 90 

*“  1  0  3  •  f  >' 

-1 02. 09  -111. 31  -119. 94 

-4.  50 

-179.88 

-116. 06 

-1 00. 52 

-102.51  -112.75  -121.36 

-5.  00 

-90. 00 

-91.91 

-96. 49 

-102.70  -114.11  -122.72 

-5.  5  ci 

-.  12 

-67. 70 

—92 . 36 

-102.80  -115.43  -124.03 

-6.  00 

-.  06 

-51.71 

-88.77 

-103.01  -116.72  -125.29 

.  50 

-.  05 

-43. 08 

-86.25 

-103.47  -118.02  -126.50 

-7.  00 

-.  04 

-39. 09 

-85. 12 

-104.31  -119.34  -127.68 

-7.50 

-.  04 

-38.24 

-85.56 

-105.62  -120.71  -128.82 

-8.  00 

-.  04 

-40. 09 

-87.65 

-107.44  -122.13  -129.92 

-8.50 

-.  05 

-45.11 

-91.37 

-109.77  —123.61  -130.99 

TABLE  3  "  TRANSFER  FUNCTION  EVALUATION  IN  THE  LAPLACE  DOMAIN 
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CONFIGURATION:  LANDS  1-C  LOW  ORDER  EQUIVALENT  L  ALF  FREE 


TRANSFER  FUNCTION  GAIN 


SIGMA! OMEGA 

.  0010 

1 . 0000 

0000 

5. 0000 

8 .  0  0  0 1 j 

1 

.  0  0  0  0 

.  0000 

.  9835 

1 . 1 397 

.  59 1  £ 

.  3735 

.  £378 

1911 

5000 

2 .  4282 

1.6158 

.  6275 

.  38£  0 

.  £399 

19££ 

- 1 .  0  0  0  0 

£4.2751 

£.1171 

.  65  03 

.  3870 

.  £41£ 

1928 

-1.50 0  0 

7 . 0375 

£. 0715 

.6545 

.  388 1 

.2415 

193  0 

-2. 0000 

£ . 6344 

1 . 6535 

.  6394 

.  385£ 

.  £408 

1 926 

-£.50 0 0 

1.6215 

1 . £7  08 

.  6092 

.  3785 

.  2392 

1918 

-3. 0000 

1. 1536 

1 . 0  0  08 

.  57  0 1 

.  3688 

•  2  o  6 1 

1905 

—3.5 0 0  0 

.  8929 

.8147 

.  5£76 

•  3566 

.  2334 

1833 

-4. 0000 

.  7273 

.  68£5 

.  46:56 

.  3428 

' .££95 

1 86? 

-4 . 5  0  0 0 

.6131 

•  5353 

.  4463 

.  328 1 

.  ££49 

1342 

-5. 0000 

.5297 

.5113 

.  41  06 

.  3 1 3 1 

.  £  1 99 

1315 

-5.5  0 0 0 

.  4661 

.  4534 

.  3785 

.  £982 

.2146 

1784 

-b. 0000 

.4161 

.  4069 

.3501 

•  233? 

.2090 

1752 

-b .  5  0  0  0 

OTtrc. 

•  '»*  ♦  j'j 

.  3689 

.  3£50 

.  £698 

.2032 

1717 

-7 .  0  0  0  0 

.  3425 

.  3373 

•  3027 

.  £566 

.1974 

1632 

-7.5  0  0 0 

.3147 

.  3 1  06 

.  £8£9 

.  £44£ 

.1916 

1646 

-3. 0000 

.£910 

.  £878 

.£653 

.  £3£6 

.  1858 

1609 

-8.5  0 0  0 

.  £706 

.  £68  0 

.  £495 

.££17 

.  1801 

1571 

TRANSFER  FUNCTION  PHAS 

E  LAG 

SIGMA?  OMEGA  . 0  0 1 0 

1 . 0  0  0  0  3 .  0  0  0  0  5.  0  0  0  0  8 .  0  0  0  0  10.  0  0  0  0 

.00  .01 

-31 . £1 

-70.57 

-84,71 

—96 . 93 

-1 03.  16 

-.50  179.05 

-44.20 

-79. 19 

-90.23 

- 1 0 0 . 46 

-106. 0 0 

-1.00  -.34 

-71.71 

-88. 69 

-95.95 

-104. 04 

- 1  08 . 86 

-1.50  -179.83  - 

1  04 .  09 

-98. 57 

-101.75 

-107. 64 

-111.74 

-£.00  -179.93  - 

128. 00 

-108. £4 

-107.53 

-111.24 

-1  14.61 

-2.50  -179.96  - 

142. 67 

-117.18 

-113. 16 

-114.79 

—117. 46 

-3.00  -179.97  - 

151.71 

-125. 12 

-118.54 

-1 18.29 

-120.29 

-3.50  -179.98  - 

157.60 

-131.97 

-123.60 

-121.71 

-123. 06 

-4.00  -179.98  - 

161. 68 

-137. 81 

-128. 30 

-125. 01 

- 1  £5 . 79 

-4.50  -179.98  - 

164. 63 

-142. 74 

- 1 32 . 62 

-128.20 

-128.44 

-5.00  -179.99  - 

1 66 . 87 

-146.92 

-136. 56 

-131.25 

—131. 02 

-5.50  -179.99  - 

168. 61 

-150.46 

-140. 14 

-134. 16 

-133. 52 

-6. 00  -179.99  - 

170.0  0 

-153.49 

-143.39 

—  1 36 . 93 

-135.94 

-6.50  -179.99  - 

171. 13 

-156. 09 

-146. 32 

-139.55 

—  1 38 . 26 

-7.00  -179.99  - 

172. 08 

-158.35 

-148.98 

-142.02 

-140. 49 

-7.50  -179.99  - 

172.88 

-160.31 

-151.38 

-144. 35 

-142.62 

-8.00  -179.99  - 

173. 56 

-162. 03 

-153.56 

-146.55 

-144.66 

-8.50  -179.99  - 

174.  14 

-163. 55 

-155.54 

-1  48.  61 

-146.61 

TABLE  4  -  TRANSFER 

function  evaluation  in  the 

LAPLACE  OOMAIN 
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LAHOS  1-C  CONFIGURATION 
^=t  rad. /sec. 
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CONFIGURATION:  LAHOS  6-2  HIGH  ORDER  SYSTEM 


TRANSFER  FUNCTION  GAIN 


SIGMA5 OMEGA 

.  0  0 1  0 

1. 0000 

.  0000 

.  9950 

1 . 3327 

— .  5  0  0  0 

.  6200 

1 . 3273 

-1.0  0  0 0 

2. 7545 

c! .  63 1 3 

-1.5 0 0  0 

.  4239 

1 . 9983 

0000 

.  0003 

a  886':' 

-2.5  0  0  0 

•  0632 

.  8 1 33 

-3. 0000 

.  4933 

1 .  0400 

-3.5  0  0  0 

1.3070 

1 . 5859 

-4. 0000 

2.9369 

2.5094 

-4.5  0  0  0 

7. 7940 

3. 3597 

-5 .  0 0 0  0-* 6 5 . 6255 

5. 1423 

-5. 5000 

1 1 . 763  0 

5.4071 

-* . 0000 

7. 0070 

4. 9886 

—6 .  5  0  0  0 

5. 5404 

4. 5611 

-7. 0000 

4.9563 

4. 3043 

-7. 5000 

4. 3063 

4.2263 

-8. 0000 

5. 0060 

4 . 3233 

-3.5  0  0 0 

5. 6 AS A 

4.6112 

3. 0000 

5. 0000 

3. 0000 

1  0 .  0  0  0 

1.2165 

.7874 

.50  05 

.3951 

1.4169 

.3541 

.5310 

.  4166 

1 . 5894 

.9183 

.5616 

.  4332 

1 . 6637 

.  9730 

.5920 

.  4593 

1 . 6576 

1 . 0324 

.  62 1 6 

.  48 1 2 

1 . 6292 

1 . 0824 

.  6503 

.  5022 

1 . 6400 

1. 1295 

a  *  l'  %  * 

.  5226 

1 . 7053 

1.175  0 

.  7  033 

.  5422 

1.8154 

1.2191 

.  7370 

.  56  08 

1 . 9483 

1.26  03 

a  r  4 Hi: 

.5781 

2. 0801 

1 . 2983 

a  766* 

.5939 

2. 1876 

1 . 3293 

.  7819 

.  6  03  0 

2 . 26  03 

1.3515 

.  7935 

.  62  01 

2.2934 

1  •  -*  6  .C'  c! 

.  8  0 1 2 

.  63  01 

2 . 3 03 1 

1 . 36 3 1 

.  8046 

.  6  379 

2.2951 

1. 3504 

•  3  036 

.  6432 

£.261 3 

1 • U247 

.  798  0 

.  6459 

2. 2031 

1 . 2859 

.  7879 

.  646  1 

TRANSFER  FUNC 
sigma; omega 

TION  RHA3 

.  00 1  0 

E  LAG 

1 .  0  0  0  0  3 .  0  0  0  0 

.  00 

.  01 

-4.  00 

-62. 66 

- .  5  0 

.  15 

1 1 . 96 

—71. 3d 

-1.0  0 

.51 

50.  18 

-82. 43 

-1.50 

.20 

1 23 . 64 

-94 . 1 8 

-2.  00 

90.19  - 

176. 28 

-103.77 

-2.50 

-.44  - 

1 26 . 43 

-110. 16 

-3.  00 

-.  14 

-93. 13 

-1  14.47 

-3.50 

-.  10 

-73. 29 

-118.36 

-4.  00 

- .  09 

—  »•"  t  «  46 

-123. 05 

-4.50 

-.  14 

-89. 03 

-129. 05 

-5.  00 

-90. 02  - 

111.79 

—1 36. 36 

-5.50 

-179.91  - 

136.29 

-144.57 

-6.  00 

-179.96  - 

153. 38 

-153.24 

-6.50 

-179.98  - 

165.31 

-162. 02 

-7.  00 

-179.99  - 

173.50 

-170.77 

-7.50 

180. 00 

179.56 

-179.52 

-3.  00 

179.99 

172.43 

171.58 

-3.50 

179.98 

163.83 

1.68.2:9 

5. 0000 
-88.25 
-93.64 
-99.33 
•105.  13 
111. 01 
116.75 
122.33 
127.95 
1 33 . 56 

139. 23 
145.17 

151.23 
157.43 
163.74 
170.  1 1 
176.50 
177.14 
It  0. 35 


8 .  o  o  o  o  i  o .  o  o  o  i; 

112. 53 

-125. 59 

•1 16.23 

-123. 75 

•120.14 

-131. 96 

124. 09 

-135.23 

1 28 .  1 2 

- 1  33 . 5 3 

1 32 . 2 1 

-141 . 85 

1 36 . 34 

-145. 19 

140. 49 

-143. 52 

1 44 . 66 

-151.34 

1 48 . 83 

-155. 12 

152.99 

-1 53. 36 

157.13 

-161.5J 

161.22 

- 1 64 . 63 

165.25 

-167. 63 

169.20 

-170.51 

1 73 . 03 

-173. 26 

1 76. 73 

-175.85 

179. 75 

-178.27 

TABLE  5  *  TRANSFER  FUNCTION  EVALUATION  IN  THE  LAPLACE -rOMAI N 
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CONFIGURATIONS  LhHOS  6-2  LOW  ORDER  EQUIVALENT  L  ALF  FIXED 


transfer  function  gain 


-•  I  SNA  •  ONEGA 

.  0  0 1  0 

1 . 0000 

3. 0000 

5. 0000 

8. 0000 

1  0 .  0  0  0 

.  0000 

.  8696 

1.4911 

1 . 2555 

.8  04  0 

.  5 1  02 

.  4  093 

-.50 0  0 

.  453 1 

2 .  1 6 1  u 

1.4131 

.  8357 

.5175 

.413  0 

-1 . 0000 

.  3952 

5 . 4836 

1 . 5526 

.  3595 

c*  I*  I-  >Zi 
m  -w'l ICO 

.4157 

—1.50 0 0 

2. 6785 

14.8261 

1 . 6260 

.  8729 

.  5259 

.4172 

-2. 0000 

3.  3  055 

4. 7201 

1 . 6-078 

.  8746 

c- 

•  C  1 

.4177 

—2.5  0  0  0 

2. 9459 

3. 0823 

1.5182 

.  8647 

.5251 

.  4 1  7  0 

-3. 0000 

2.4183 

2. 3457 

1 . 3962 

.  8443 

.5213 

.4152 

-3. 5000 

1.9321 

1 . 3958 

1 . 2697 

.8174 

.  5154 

.4123 

-4. 0000 

1 . 6532 

1.5350 

1.1519 

.  735  0 

.  5077 

.4084 

-4 .  5  0  0  0 

1 . 4067 

1 . 3566 

1 . 0468 

.  7499 

.  498  3 

.  4036 

-5. 0000 

1.2137 

1 . 1323 

.  9546 

.7139 

.  4377 

.  393  0 

-5 . 50  0  0 

1 . 0722 

1 . 0454 

.3743 

.  678 1 

.  47 62 

.  39 1 3 

-A . 0  0  0  0 

.  9555 

.  9354 

.  3044 

.  6435 

.  46  39 

.  3849 

-6.5  0 0 0 

.  36  0  7 

.  3454 

.74  33 

.6105 

.4512 

.  3777 

-7. 0000 

.7324 

.  77  05 

.  6397 

.  579  3 

.  4332 

.  3  7  0  0 

-7. 5000 

•  (■'  1 6  - 

.  7074 

.  6425 

.  55  0 1 

.  4252 

.  362 1 

-3. 0000 

.  661  0 

.  6535 

.6007 

.  5229 

.4122 

.  354  0 

-8 . 5  0  0  0 

.6131 

.6070 

.  5636 

.  4977 

.  3994 

.  3459 

TRANSFER  FUNCTION  PHA 

3E  LAG 

i I  UNA  ? OMEuA  . 0  0 1 0 

1 . 0  0  0  0 

3. 0000  5. 0000 

3. 0000 

10.  0  0  0 

.  00  . 02 

-3.53 

-64.09  -90.49  - 

114. 00 

—  1 26 . 58 

-.50  .21 

1 1 . 34 

-72.53  -96.14  - 

117. 58 

-129. 44 

-1.00  179.76 

34 . 84 

—83.42  —1  Oil.  22  — 

121.25 

-132. 35 

-1.50  179.94 

—  1 82  •  46 

-96.04  -108.57  - 

124.98 

-135.29 

-2. 00  -130. 00 

-157.67 

-108.77  -115.00  - 

1  c8 .  i*‘  3 

-138.25 

-2.50  -179.98 

—  1 5 1 '  .  c  0 

-120.17  -121.33  - 

132.48 

-141.20 

-3.00  -179.98 

-159.73 

-129.71  -127.38  - 

1 36 .18 

-144.  13 

-3.50  -179.98 

-162. 67 

-137.49  -133.05  - 

1 39 . 80 

-147. 02 

-4.00  -179.99 

-165. 35 

-143.84  -138.27  - 

143.31 

-149. 86 

-4.50  -179.99 

—  1 6 1  .  6  c 

-149.03  -143.02  - 

146.71 

-152.64 

-5.00  -179.99 

-169.51 

-153.45  -147.32  - 

1 49 . 96 

-155. 35 

-5.50  -179.99 

-171.09 

-157.13  -151.20  - 

1 53 . 06 

-157.97 

-6.00  -179.99 

-172.40 

—160.27  —154.69  — 

156. 00 

-160.51 

-6.50  -179.99 

-173.50 

—162.97  —157.83  - 

158.79 

-162.95 

-7.00  -179.99 

-174. 44 

-165.32  -160.66  - 

161.41 

-165.29 

-7.50  -180.00 

-175.24 

—167.36  -163.22  - 

163. 88 

-167.53 

-8. 00  -180. 00 

-175.94 

-169.16  -165.53  - 

166.21 

-169. 68 

-8.50  -130.00 

-176.54 

-170.74  -167.62  - 

1 68 . 39 

-171 . 73 

TABLE  6  -  TRANSFER  FUNCTION  EVALUATION  IN  THE  LAPLACE  DOMAIN 
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Figure  3  compares  the  forced  response  of  the  LAHOS  1-C  configuration 
with  its  first  over  second-order  with  time  delay  equivalent. 

At  c  =  -1.5: 

Ma9H0S  =  '-°97  eH0S  =  ■'97'' 

Ma9L0S  =  2.07!  -LOS  = 

which  constitute  substantial  magnitude  and  phase  errors. 

Figure  18  represents  the  total  time  response  comparison  of  these  two 
systems  for  the  input, 

x(t)  =  Aect  sin  jt,  3  =-1.5,  -  =  1.0 
and  shows  that  of  the  high-order  and  low-order  systems  to  be  highly  similar. 

Recall  that  in  setting  the  similarity  conditions  on  the  total  time 
response,  G  and  3  were  considered  variables.  If  they  were  set  equal  to 
the  high-order  system  values,  total  variables  and  therefore  similarity 
conditions  would  be  lost. 

Therefore,  it  is  clear  that  similarity  of  forced  response  is  not  a 
general  principle  that  holds  over  the  S  plane  nor  in  fact  should  this 
condition  be  enforced  in  determining  equivalents.  Only  the  total  time 
response  can  be  the  comparison  standard  over  the  S  plane.  It  will  be 
shown  later  that  importance  of  the  forced  response  along  the  ju  axis 
lies  in  the  fact  that  it  also  represents  the  total  long-term  or  steady 
state  response  to  the  general  input  at  those  points  in  the  S  plane. 
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TIME  HISTORY  COMPARISONS 


Because  the  low-order  system  equivalents  were  obtained  purely  by 
matching  the  system  forced  response  along  the  line  r=  0,  it  is  likely 
that  the  six  time  response  similarity  conditions  are  not  met  at  all 
points  within  the  S  plane  region  of  interest,  even  though  good  frequency 
response  matches  have  been  obtained.  To  determine  the  extent  of  the 
mismatching,  time  history  responses  were  calculated  for  input  points 
across  the  defined  S  plane  region.  Figure  5  depicts  both  the  region 
correspond i ng  to  a  maximum  frequency  of  10  radians/second,  and  a  time 
to  init:ai  input  peak  of  tQ  =  .157  seconds  and  the  distribution  of 
input  test  signals  across  that  region.  Table  7  relates  the  frequency 
damping  and  amplitude  of  each  input  to  a  letter  designation  for  easy 
reference . 


CASE  NUMBER 

J 

0 

A 

A 

0.0 

1 

.0 

1.0 

8 

"1.5 

4.355 

C 

-4.0 

10.985 

0 

-6.0 

16.385 

E 

0.0 

5 

.0 

1  .0 

F 

-1  .0 

i 

1  .342 

G 

-3.0 

2.164 

H 

-5.0 

3.102 

1 

0.0 

s’ 

.0 

1  .0 

■*> 

-0.5 

1.10 

K 

-1.5 

1.319 

L 

-2.5 

1 

1.557 

x(s)  =— 

(a  +  juj)  J  Is 

- 

(a  -  ju)] 

e-a/w  tan" 1 (~£) 

A  - - - - 

sin  { tan" 1  (— )  } 

a  ' 

TABLE  7:  TEST  INPUT  CASES 
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Conditions  of  constant  unity  magnitude  on  the  initial  peak  of  the 
input  signal  were  imposed  in  order  to  prevent  the  input  signal  from  dimin¬ 
ishing  in  magnitude  as  damping  was  increased.  Time  for  the  initial  peak 
has  been  shown  to  be: 


1  f  the 


input  value  at  this 


time  is 


unity: 


x(tQ)  =  Ae°*’°  sin  ^tD  =  1 


A  = 


-o/ui  tan-1  (— ) 
e _ o 

sin  {tan*1  (-y) } 


The  time  responses  were  calculated  using  a  digital  computer  program  which 
employs  the  Heaviside  expansion  to  calculate  inverse  Laplace  transforms 


y(t)  =<£~1{F(s)  x(s)| 

where  F(s)  is  the  transfer  function  and  x(s)  =  -r-^—^r 

\~a  ,jj) 

Since  the  author  was  not  responsible  for  this  program,  its  results 
were  verified  by  a  hand  solution  of  the  time  response  to  input  case  B 
for  the  LAHOS  1 -C  equivalent.  This  verification  appears  in  Appendix  B. 

The  LAHOS  1-4  conf i gura t i on  was  run  for  the  entire  range  of  test  sig¬ 
nals.  The  results  appear  in  Figures  6  through  17-  Because  the  total  out¬ 
put  magnitude  decreases  with  both  uj  and  a,  the  absolute  importance  of  the 
response  differences  diminishes  at  higher  frequencies.  It  was  decided  to 
examine  the  LAHOS  1-C,  and  LAHOS  6-2  configurations  only  at  the  lowest 
frequency  cases  where  the  differences  between  high  and  low  order  system 
responses  are  most  noticeable.  These  appear  in  Figures  18  through  23. 

Figures  6,  10,  and  14  show  the  LAHOS  1-4  configuration  response  to 
the  undamped  sine  wave  inputs.  The  responses  show  excellent  agreement 
through  the  initial  peak  which  progresses  toward  small  errors  in  the 
magnitude  and  phaseing  of  the  stady  state  response  as  indicated  by  the 
small  differences  in  the  frequency  responses  of  reference  1.  As  we  pro¬ 
ceed  into  the  S  plane  at  each  frequency,  slight  decreases  in  the  accuracy 
of  the  initial  response  peak  may  be  detected  along  with  a  quite  noticeable 
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LAHOS  J-4  CONFIGURATION 


FIGURE  6  -  Pitch  Rate  Response 


LAHOS  1-4  CONFIGURATION 


INPUT  3 
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LAHOS  1-4  CONFIGURATION 


unc 

LOS - 

— 

£ 

\  cl 

«H 

H - 

0  2  4-  6  3  10 


t  ^  sec. 


FIGURE  9  ■  Pitch  Rate  Response 
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LAHOS  1-4  CONFIGURATION 


t  ^  sec 


Figure  11  -  Pitch  Rate  Response 


MADC - 7923 1 -60 


LAHOS  1-4  CONF I GURATICN 


0  12  3  4  5 


t  ^  sec 


Figure  13  -  Pitch  Rate  Response 
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LAHOS  1-A  CONFIGURATION 


Figure  14  -  Pitch  Rate  Response 
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LAHOS  1-4  COMF1 GURAT i ON 


t  ^  sec . 


FIGURE  15  "  Pitch  Rate  Resoonse 
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INPUT  K 


,  i 

o  I 


2  3 


t  ^  sec. 


FIGURE  16  -  Pitch  Rate  Response 


INPUT  3 


FIGURE  1 8  -  Pitch  Rate  Response 
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LAHOS  6-2  CONFIGURATION 


0  2  4  6  3  10 
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LAHOS  6-2  CONFIGURATION 


0  2  '■*  6  3  10 

t  ^  sec 


Figure  22  -  Pitch  Rate  Response 
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FIGURE  23  '  PITCH  RATE  RESPONSE 
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difference  in  the  subsidence  times  of  the  high  order  system  (tH)  and 
low  order  system  ( t l ) .  The  relative  trend  persists  at  high  frequencies 
but  is  less  noticeable  due  to  the  decreased  magnitude  of  the  output. 

Figures  18  to  20  show  the  LAHOS  1-C  high  and  low-order  system  re¬ 
sponses  to  the  low  frequency  well-damped  inputs.  The  higher  order 
equivalent  shows  a  much  better  overall  match  of  the  time  response  by 
improving  the  subsidence  time  match  while  incurring  slightly  greater 
initial  peak  errors.  Similar  statements  hold  for  the  time  responses  of 
LAHOS  6-2  in  Figures  21  to  23. 

While  it  may  be  doubtful  that  these  differences  (for  LAHOS  1-6, 
and  LAHOS  6-2)  would  be  significant  to  the  pilot  observing  them,  they 
do  show  that  the  equivalent's  ability  to  represent  the  high-order  sys¬ 
tem  varies  with  the  input  signal  and  is  not  optimized  for  any  input  or 
S  plane  point.  The  LAHOS  1-4  configuration,  while  not  using  the  recom¬ 
mended  form  for  the  equivalent,  demonstrates  that  j'jj  Bode  matching  alone 
does  not  guarantee  adequate  similarity  of  response.  Further,  the  theory 
prediction  that  this  form  could  not  satisfy  all  the  similarity  conditions 
is  confirmed. 
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As  the  damping  of  the  input  signal  is  decreased  to  zero,  the  long¬ 
term  output  does  not  attenuate  but  approaches  the  forced  response.  Sim¬ 
ilarity  of  the  forced  response  to  undamped  inputs  is  the  same  condition 
as  matching  the  subsidence  time  for  well-damped  inputs.  They  both  match 
the  long-term  or  steady-state  response. 

The  six  non-linear  equations  which  enforce  the  total  response  sim¬ 
ilarity  conditions  are  dependent  on  the  value  of  S;  that  is,  the  input 
chosen.  In  general  the  equivalent  systems  parameters  will  vary  over 
the  Laplace  domain.  The  exception  is  the  case  where  the  high-order 
system  and  its  equivalent  are  of  the  same  form. 

In  the  case  where  the  equivalent  system  parameter  variation  over 
the  region  of  interest  is  small,  they  might  be  considered  constants. 

The  parameters  could  equally  well  be  determined  by  satisfying  a  single 
similarity  condition  for  several  different  inputs  as  by  satisfying  sev¬ 
eral  similarity  conditions  for  a  single  input.  To  within  acceptable  ac¬ 
curacy  the  same  parameters  would  be  obtained.  Locating  equivalent  sys¬ 
tem  parameters  by  matching  jw  Bodes  in  effect  assumes  they  are  approx¬ 
imately  constant  over  the  S  plane  and  imposes  the  similarity  of  the  long 
term  response  along  the  line  a  =  0.  The  initial  similarity  conditions 
are  always  matched  in  the  determination  of  the  free  response  mode  coeffi¬ 
cients.  If  a  very  good  frequency  response  match  is  obtainable,  the  as¬ 
sumption  of  invariance  is  good  at  least  along  the  jw  axis,  and  adequate 
matching  of  other  similarity  conditions  (with  sinusoidal  inputs)  is  implied. 
As  shown  by  the  LAHOS  1-4  time  histories,  similarity  of  the  j m  Bode  alone 
does  not  guarantee  time  response  similarity  elsewhere. 

One  method  of  obtaining  time  response  matched  equivalents  has  been 
suggested  here.  Others  might  be  possible.  The  important  point  is  that  the 
equivalent  system  parameters  vary  with  both  a  and  u  and  that  j*  Bode 
matching  techniques  depend  on  that  variation  being  negligible.  This  im¬ 
plication  remains  regardless  of  the  method  used  to  match  the  time  response. 
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Failure  to  obtain  an  acceptable  8ode  match  is  a  direct  indication  that 
the  equivalent  system  parameters  vary  significantly.  The  requirements  of 
MIL-8785B  allow  these  parameters  to  vary  within  certain  tolerances.  The 
fact  that  they  vary  does  not  mean  that  they  exceed  the  specified  toleran¬ 
ces.  Failure  to  obtain  a  good  jw  Bode  match  is  therefore  an  insufficient 
reason  for  discarding  the  system. 

It  is  too  restrictive  to  require  that  for  a  system  to  be  judged  ac¬ 
ceptable  that  a  single  low-order  equivalent  be  capable  of  representing  it 
over  the  entire  S  plane  region.  Rather,  it  is  necessary  only  that  the 
high-order  system  behave  like  some  acceptable  lower-order  system. 

The  success  of  jm  Bode  matched  equivalents  here  and  in  the  references 
indicates  that  in  many  cases  invariance  of  the  equivalent  system  parameters 
is  a  good  assumption.  But  uncertainty  for  each  individual  case  will  always 
exist  until  the  assumption  is  checked.  Further,  time  delays  are  required 
in  many  cases  to  enforce  the  invariance  assumption.  As  discussed  on  page 
17,  this  alters  the  parameters  and  decreases  the  similarity  of  the  time 
response.  It  is  uncertain  that  a  low-order  system  with  a  small  time  delay 
will  always  be  evaluated  in  the  same  manner  as  an  identical  system  without 
any  time  delay.  Yet  that  assumption  is  made  in  applying  the  current  equi¬ 
valent  system  parameters  to  the  requirements  of  MIL-8785B.  If  the  time 
delay  is  sufficiently  small  this  is  not  a  serious  objection,  but  no  such 
difficulties  are  encountered  with  time  response  matching. 

Section  3.2.2. 1  of  MIL-8785B  describes  tolerances  for  Csp  and  ^nSp 
as  a  function  of  Nz/a  (La  or  I/T92)  Given  a  low  frequency  magnitude  of 
0  db  and  low  frequency  phase  angle  of  0°,  these  tolerances  define  a  magni¬ 
tude  and  phase  envelope  in  the  frequency  domain: 


Envelope 


*  (  Not  plots  of  actual  envelopes  ) 
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If  a  high-order  system  has  a  frequency  response  within  that  envelope,  then 
there  is  some  acceptable  low-order  system  which  has  the  same  steady  state 
response.  If  the  high-order  system  response  violates  the  envelope  it  may 
be  discarded.  If  it  possesses  a  good  lower  order  system  match,  the  equi¬ 
valent  system  parameters  are  approximately  constant  for  undamped  inputs 
within  the  frequency  range.  Significant  variations  may  still  exist  in  the 
a  direction.  If  the  high-order  system  response  lies  within  the  envelope 
but  does  not  possess  an  acceptable  frequency  response  matched  equivalent 
the  only  way  to  determine  the  equivalent  system  parameters  and  therefore 
its  acceptab i 1 i ty  is  by  total  time  response  matching. 

Current  Bode  matching  techniques  are  a  good  test  of  the  high-order 
system  only  if  the  implied  assumption  of  parameter  invariance  holds  over 
the  S  plane.  The  assumption  works  well  in  many  cases  as  long  as  time  delays 
are  included  in  the  equivalent  system  form,  but  statistical  data  can  only 
establish  the  probability  of  its  holding  in  any  given  application. 


IflECiiilw;  PjiGE  Bx^diK  -  i.cr  F1^±D 
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IMPLICATIONS  TO  PILOT  COMPENSATION  CRITERIA 


Reference  3  developes  a  flying  qualities  criteria  based  on  pilot  com¬ 
pensation  of  the  form, 


os  (TjS  +  1) 

6c  "  P  Tf2S  +  2) 


e 


-  .  3S 


The  basic  philosophy  is  that  pilots  desire  a  certain  'optimum"  response 
and  are  willing  to  provide  limited  amounts  of  compensation  to  achieve  that 
response  before  down-rating  the  system.  The  final  pilot  rating  is  then 
dependent  on  how  much  compensation  must  be  provided  and  the  nearness  of 
the  best  result  to  the  optimum. 

If  it  is  assumed  the  pilot  wants  direct  pitch  control  with  no  phase 
lags  the  implication  is  that  at  low  frequency: 

Mag  (6/6c)  =  1 
Mag  (e/9c) Qb  =  0 
-4“  (a/3c)  =  0 

Deviations  from  this  (droop)  should  be  minimized. 

Furthermore,  if  oscillatory  overshoot  is  to  be  minimized,  a  minimum 
zero-over-second-order  form  for  the  optimum  must  be  considered  in  order 
for  the  system  to  have  an  oscillaroty  free  response  mode.  Maintaining  a 
certain  bandwidth  maintains  the  desired  response  up  to  a  certain  speed  of 
input.  For  a  well-damped  (5  =  .707)  0/2nd  system  the  -3  db  bandwidth  oc¬ 
curs  at  the  break  frequency,  ,  where  the  phase  angle  passes  through  -n/2. 
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The  compensation  parameters  Kp,  T-,  T ; ,  are  determined  by  providing 
phase  compensation  to  hold  t  =  -*/2  at  the  specified  bandwidth,  while 
matching  the  low  frequency  magnitude  and  minimizing  droop.  The  amount 
of  phase  compensation  is  then  a  measure  of  how  hard  the  pilot  must  work 
to  obtain  the  optimum  and  the  resonance  is  a  measure  of  how  close  to  the 
optimum  he  is  able  to  get.  Since  the  object  is  to  come  as  close  as  pos¬ 
sible  to  the  simple  "optimum"  Bode,  this  method  is  providing  the  same  re¬ 
sults  as  matching  the  system  Bode  to  obtain  an  equivalent.  In  this  case 
we  are  attempting  to  match  the  aircraft  and  flight  control  system  dynamics 
plus  pilot  compensation  with  the  optimum  zero  over  second  order  system. 

The  main  difference  is  that  we  are  attempting  to  match  the  low-order  sys¬ 
tem  by  varying  root  locations  in  the  high-order  system  rather  than  vice- 
versa  . 

If  we  view  this  criteria  in  terms  of  the  high-order  system  Bode,  we 
see  that  by  allowing  the  pilot  to  provide  certain  magnitude  and  phase  com¬ 
pensation,  the  high-order  system  magnitude  and  phase  may  be  allowed  to  vary 
from  the  optimum  by  a  like  amount.  This  is  the  same  as  allowing  the  high- 
order  system  Bode  to  lie  within  a  certain  envelope  similar  to  that  already 
discussed.  This  technique  then  is  not  as  restrictive  as  current  j  Bode 
matching  techniques  because  it  allows  magnitude  and  phase  to  vary  within 
that  envelope. 

The  pilot  compensation  approach,  however,  is  sensitive  to  the  band¬ 
width  chosen.  Like  the  current  equivalent  systems  methods  it  does  not  con¬ 
sider  what  differences  in  total  response  to  highly  damped  inputs  might  oc¬ 
cur.  Thus  the  pilot  compensation  approach  is,  in  principle,  doing  much  the 
same  thing  as  current  equivalent  systems  approaches  and  suffers  from  some 
of  the  same  deficiencies. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


Current  methods  of  applying  equivalent  systems  as  well  as  proposed  pilot 
compensation  criteria  are  based  on  sound  concepts.  Both  have  been  shown  to 
have  highly  similar  implications  and  some  common  limitations.  Current  equi¬ 
valent  systems  methods  are  limited  by  the  implied  assumption  of  invariance 
of  the  equivalent  systems  parameters.  Since  they  deal  only  with  Bodes 
neither  method  considers  total  system  response  to  aperiodic  pilot  inputs. 
Provided  the  recommended  equivalent  system  form  which  includes  artificial 
time  delays  is  employed,  the  invariance  assumption  is  valid  in  many  cases. 

But  because  it  is  not  assured,  straightforward  application  of  either  method 
may  not  always  be  reliable.  Determination  of  equivalent  systems  para¬ 
meters  based  on  time  response  matching  not  only  eliminates  the  uncertainty, 
but  is  directly  implied  both  by  the  equivalent  system  concept  and  by  the 
idea  of  an  optimum  response  basic  to  pilot  compensation  methods.  As  such 
it  represents  the  most  logical  extension  of  current  specifications. 

The  analysis  suggests  that  equivalent  systems  parameters  be  considered 
variables  over  the  region  of  the  Laplace  domain  within  which  human  pilots 
may  operate,  that  their  values  be  determined  by  similarity  of  total  time 
responses  to  inputs  within  that  range;  and  that  the  acceptability  of  the 
system  be  judged  by  comparing  the  variations  with  the  current  MIL- F-8785B 
requ i rements . 

It  is  recognized  that  application  of  this  criteria  presents  some 
problems.  Programming  simulators  for  direct  comparison  of  high-order 
systems  with  their  equivalents  is  difficult  because  of  the  variance  of 
equivalent  systems  parameters  with  input.  No  design  methods  currently 
exist  to  meet  the  criteria.  It  is  unclear  what  judgement  should  be  made 
about  a  high-order  system  if  it  is  acceptable  over  only  a  portion  of  the 
S  plane  region.  Additional  verification  of  the  implications  are  necessary 
before  high  levels  of  confidence  in  the  method  are  achieved.  In  particular 
some  of  the  anomalies  which  have  occured  with  present  methods  should  be 
examined  under  the  suggested  criteria. 

It  is  believed  that  most  of  these  problems  may  be  either  solved  or 
avoided  as  further  attempts  to  apply  the  methods  are  made. 
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APPENDIX  A 

APPLICABLE  PROPERTIES  OF  DIFFERENTIAL  EQUATIONS 


A-l 


fHEOsaiwo  Pilcar  _ 


«cr 


NADC-79231-60 


All  linear  control  theory  considers  solutions  to  a  linear  ordinary 
differential  equation  with  constant  coefficients: 


n 

Z 

i  =0 


djy (t;  _ 
dt1 


dix(t) 

dtJ 


(A-  1 ) 


where  the  output  function,  y(t),  is  initially  zero  and  for  times  greater 
than  zero  is  a  consequence  of  the  input  function,  x(t).  For  any  given 
input  function  there  may  be  many  solutions  for  the  output  function,  the 
complete  solution  being  the  sum  of  all  such  solutions.  It  is  a  direct 
consequence  of  the  system's  linearity  that  any  sum  of  solutions  for  the 
output  is  also  a  solution  and  that  the  complete  response  to  any  linear 
combination  of  inputs  is  the  sum  of  the  complete  responses  to  each  in- 


d  i  v  i  d  ua  1 

input. 

IF: 

y  1  ( t) 

and  y 2 ( t )  are  solutions, 

THEN 

:  yi(t) 

+  y2(t)  is  also  a  solution. 

and 

IF: 

Xi  (t) 

produces  y  1  (t) , 

and : 

x2  (t) 

produces  y2  (t)  , 

THEN 

X  1  ( t ) 

+  x2(t)  produces  y2(t)  +  y2(t) 

S i nee  zero  may  be  added  to  either  side  of  equation  (A-l)  without 
changing  the  equality,  any  functions  which  make: 


n 


Z  Bj 

i  =0 


d'y  (t) 
dt1 


=  0 


(A-2) 


m 


j*0 


A .  dJ.*.U).  =  0 

J  dtJ  U 


(A-3) 


are  always  solutions.  The  functions  for  y(t)  implied  above  are  called  the 
homogeneous  or  free  response  solutions  since  they  do  not  depend  directly  on 
the  input  function  (s).  All  other  solutions  depend  explicitly  on  the  input 
function  (s)  and  are  termed  the  forced  response.  The  total  response  is  the 


I 


A-3 
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sum  of  both:  y (t)  =  y (t) forced  +  y(t)free_ 

Equations  of  the  form  (A-1)  are  easily  solved  using  the  Laplace 
transform  which  changes  the  operation  of  differentiation  into  multipli¬ 
cation  by  the  Laplace  variable,  S. 


/[^]  -  »'»<■ 


changes  (A-1)  to: 


n  m 

y(s)  E  BjS»  =  x ( s )  E  A-SJ 
i =0  j  =0 

which  when  factored  and  rearranged  yields  the  transfer  function: 


IlL  =  i=0  tS  -  (a;  +  j  b ; )  _/ 


{S  -  (a ;  +  ibj)} 


=  i=£ 
xtsT  Q 


The  above  function  in  the  complex  variable  S  *  a  +  jw  may  be  evaluated  as 
a  single  complex  number  at  any  point  in  the  Laplace  domain: 

44  =  GeJ* 

x(s) 

Choosing  a  value  for  s  is  equivalent  to  setting  the  operation  of  differen¬ 
tiation  equal  to  multiplication  by  a  complex  constant. 


IF: 

s  =  a  +  jui 

and  : 

»<•> 

THEN : 

/pfr1]-  <»  *  *<s> 

i nverse 

Laplace  transforming: 

~4~'  *  + 

f  ^44  =  (o  +  jw)  /  dt  +  C 

J  x(t) 
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In  x(t)  =  (a  +  j*) t  +  C 
x(t)  =  Ae(c  +  j-)c  A  =  eC 


The  transfer  function  equation  becomes: 

y{s)  =  GeJ'3  x ( s ) 
inverse  Laplace  transforming 

y  (t)  =  GeJ “  x ( t ) 


or 


y(t)  =  GAe0t  eJ'^t  +  9) 


For  x(t)  and  y(t)  to  be  real  functions  of  time,  a  and  w  must  appear  in 
conjugate  pa i rs : 

x(t)  =  Ae3t  e-i'ajt 
or 


x(t)  =  Aeat  (cos  jjt  +  j  sin  *t) 


A  pair  of  inputs  g  i  ves : 


x2(t)  +  x2  (t)  =  A1e°lt(cos  Jit  +  j  sin  jit)  +  A2e"2L  (cos  1  +  j  sin  -u2t) 
For  the  imaginary  part  to  vanish: 

A1eait  sin  Ult  +  A2e''zt  sin  a)2t  =  0 


which  directly  implies: 


A,  =  A 


2 1  a. 


and 


Xl 


(t)  +  x2(t)  =  Aeat  (eJj)t  +  e'jujt) 


ha  I  f  of  wh i ch  is: 


*i(t)-  .  x  2  ( t )  =  .at 
2  +  2  r’©  cos  cut 

Realizing  that  the  arbitrary  constant  of  integration  might  be  complex  makes 
A  in  general  a  complex  coefficient,  allowing  phasing  of  the  input. 


A 
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x(t)  =  Ae'" c  sin  Ut  +  o) 
y(t)  =  GAe' C  sin  (~t  +  5  +  i) 

represent  the  general  form  of  the  input  and  that  part  of  the  output  directly 
dependent  on  the  Input,  the  forced  response. 

If  equation  (A-2)  is  Laplace  transformed  it  is  identical  to  setting 
the  denominator  of  the  transfer  function  ( A- 5 )  to  zero.  The  roots  of  the 
character i st i c  equation  then  represent  the  free  response  solutions, 

n 

y(t)  =  I  B;er'  where  r  j  =  a  j  +  j  b  j 

i  =1 

Since  the  free  response  is  also  a  real  function  of  time  the  same  argu¬ 
ments  about  conjugate  roots  and  complex  coefficients  hold  as  for  the  for¬ 
ced  response.  The  completely  general  form  of  the  complete  response,  al¬ 
lowing  for  simultaneous  input  functions,  is: 

y(t)  =  lGAeat  sin{a)t  +  (d+S)}  +  EBeat  sin(bt  +  u) 

\ _  /  \  / 

v  - v - 

Forced  Response  Free  Response 

Since  the  free  response  represents  all  solutions  to  the  homogeneous 
and  characteristic  equations  it  is  complete.  Some  question  might  remain 
as  to  the  completeness  of  the  forced  response.  To  examine  this  question 
we  may  seek  additional  solutions  for  a  simple  case  using  the  known  forced 
response  solution. 

Consider  the  zero  over  first  order  system  described  by  the  differen¬ 
tial  equation: 

+  ay  (t)  =  x  ( t ) 

whose  transfer  function  is: 

y (s)  3  1 

x(s)  s  +  a 

Choosing  the  input  x(t)  *  e  at  is  equivalent  to  setting  s  =  -a.  The 
known  forced  solution  is  of  the  form  y(t)  *  Ge”0t. 


A-6 
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With  this  function  of  the  output  the  differential  equation  may  be  written 
as : 

-oGe“J<:  +  aGe*at  =  e  c  or  G  =  -7 — - 7 

which  is  the  magnitude  of  the  transfer  function  at  s  =  -1.  Additional 
solutions  may  be  found  by  multiplying  the  known  solution  by  another 
function  of  time,  h(t),  or  by  adding  an  additional  function  of  time,  g(t) 

y(t)  =  h ( t )  Ge'at  +  g(t) 

where  neither  h(t)  nor  g(t)  equal  Ge-a-. 

Returning  again  to  the  differential  equation, 

y 1 ( t)  =  h 1 ( t )  Geat  +  g‘(t)  -a  G  h(t)e"^t 

and 

h‘(t)  Ge_c7t  +  (a-a)G  h(t)e  ot  +  { g  *  ( t )  +  a  g(t)}  =  e~Jt 
equating  the  coefficients  of  e~0t  gives  the  homogeneous  equation: 

g1 (t)  +  a  g(t)  =0 

and 

G  {h1 (t)  +  (a-o)  h (t) >  =  1 
Making  the  substitution: 

u(t)  =  h (t)  -  1 
u' (t)  =  h' (t) 

the  second  equation  becomes: 

u *  ( t )  +  (a-o)  u ( t )  =0  -jl  =  (a-a). 

which  is  the  same  form  as  the  homogeneous  equation  and  has  the  solution: 

u(t)  =  Ce-(°'a)t 
h(t)  =  1  +  Ce'(a'a)t 

The  total  solution  is: 

y (t)  =  Ge'at  {1  +  Ce'^~a+a)t}  +  Be'at 

°r  y ( t)  *  Ge-crt  +  {GC  +  B)  e"at 


A-7 
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Applying  the  condition  that  the  output  is  initially  zero: 
y  (0)  =  G  +  (GC  +  B)  =  0 
or 

y(t)  =  {e'Jt  - 

and  both  h(t)  and  g(t)  have  collapsed  back  into  the  free  response  solution. 
With  some  additional  complexity  the  same  principle  holds  for  a  zero  over 
second  order  system.  Although  more  rigorous  proof  will  not  be  attempted 
here  it  may  be  taken  that  equation  (A-6)  represents  the  complete  solution 
except  at  the  poles  of  the  transfer  function. 

In  the  expression  above  as  c  approaches  a  the  output  becomes  undefined. 


I  im 
a-*a 


—  (e”at 
a-a 


In  this  case  we  may  use  the  free  reponse  solution  to  find  the  complete 
solution  because  it  is  always  part  of  the  complete  solution. 


V^free  =  8e  ‘ 

Additional  solutions  are: 

y (t)  =  h ( t )  Be  at  +  g(t) 

Replacing  y(t)  in  the  differential  equation  with  this  expression: 

h' (t)  Be  at  -  a  h(t)  Be  at  +  g'(t)  +  a  h(t)  Be  at  +  a  g ( t )  =  Ae  at 

which  yields  again  the  homogeneous  equation:  g'(t)  +  a  g(t)  =  0,  and 

h1 (t)  Be'at  =  Ae'at 
h1  (t)  B/A  =  1 
h ( t )  *  (A/B) t  +  c 

Then: 

y ( t )  *  (At  +  BC)e  at  +  Ae  at 
y (t)  -  (At  +  (A+BC)  }  e*at 


A-8 


For  the  initial  condition: 

y (0)  =  (A+SC)  =  0 

the  complete  solution  at  a  =  -a  is: 
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SUMMARY 


The  properties  of  linear  systems  used  in  this  analysis  are  then.- 

1.  The  total  time  response  of  the  system  is  the  sum  of  the  free  response 
and  forced  response. 

2.  The  transfer  function  magnitude  and  phase  angle  at  any  point  in  the 
Laplace  domain  represents  the  complete  forced  response: 

y(t)  =  GAev  t  sin  {oot  +  p  +  d) 
to  the  i nput  s i gna 1 : 

x(t)  =  Aeat  sin  (^t  +  P) 

3.  The  free  response  is  determined  by  the  transfer  function  pole  locations: 

y(t)  =  IBneant  sin  (bnt  +  Pn) 

and  the  complex  coefficients,  Bne|:^n  are  determined  by  satisfying 
initial  conditions  (usually  zero)  on  the  total  time  response  and  its 
f i rst  n  -  1  der i vat i ves . 

4.  The  total  time  response  to  the  sum  of  any  such  inputs  is  the  sum  of  the 
total  response  to  each  input  separately. 

5.  The  only  exception  to  the  general  output  form  is  at  the  pole  locations 
where: 

y(t)  =  Ate  a 


A- 10 
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APPENDIX  B 


CHECK  OF  INVERSE  LAPLACE  TRANSFORM  PROGRAM 


I 


B-l 


■  ERECiQlNG  PjlGE  BlMUL  -  hC  f  n  . ..  i.  n 
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The  LAHOS  1-C  equivalent  system  from  page  20  is: 


y(s)  _  1-9257  (S  +  .38)  e 
x (s)  ~  (s  +  1  . 14)  (s  +  .641 


■•037s 


For  input  case  8,  page  33: 

x(t)  =  4.355  e'1 -5t  sin  (t) 

Table  5  shows:  G  =  2.071,  5  =  -lOA0,  1 . 8 1 5  radians,  and  the  output  is: 

y(t)  =  GAejt  sin  (wt  +6')  +  Bie3^  +  B2ea:t 
To  determine  the  free  response  mode  coefficients  set, 
y  (0 )  =  GA  sin  a  +  Bj.  +  B2  =  0 
y'(0)  =  GA  (ucos6‘  +  a  sin  6')  +  B ,  a :  +  B  2a  ,  =  0 
which  may  be  written: 

GA  Lcose'  +  (a~ a  )  sin  6')  +  Bz  (a2  -  aj  =0 
GA  (ujcose '  +  (a-a2)  sin  31}  +  BL  (a>  -  a2)  *  0 


—  =  vicos»‘  +  (g-a  2)  sin  e' 

GA  ( a  ^  -  a  2  j" 

li.  a  —  WCOS91  +  (c~ai)  sin  3' 

GA  ( a2  - a2 ) 

The  above  phase  angle  includes  an  increment  due  to  the  time  delay  of: 
23  *  -jT *  .037  radians 

SO  3'  =  9  -29  ,  9'  =  T . 778  radians 

with  v=*  1,  a  *  -1 .5,  a1  =  -l.t4,  a2  =  -.6487, 

«  *  '-2717 

-  -.2935 

GA  =  9.019 
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and  Che  output  function  is: 

y(t')  =  9-019  (e’1‘5t'  sin  (t 1  -  1.778)  +  1.2717  e‘ 1 ' ' 4t- . 2935e' -6487t ' 

by  inserting  values  of  t'  =  t  -  T  =  t  -.037  sec.  into  this  equation  values 
for  y(t)  may  be  calculated.  This  was  done  and  compared  to  the  values  com¬ 
puted  by  the  program  in  Table  B-l  and  were  found  to  agree  within  the  error 
shown . 

TABLE  B-l 


time  (sec.) 

Program 

Equat i on 

Error 

1 

.8677 

.3608 

.0069 

2 

•  5^67 

.5701 

•  0239 

3 

.0900 

.1022 

.0122 

4 

-.0606 

-.0579 

.0027 

5 

-.0658 

-.0660 

.0002 

The  worst  error  is: 

e  =  .024, 

max 

and  the  average 

error 

4 
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